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>A  three-dimensional  computer  code  has  been  developed  for  the  solid-fuel 
ramjet  combustion  process.  The  code  includes  finite-rate  chemical  kinetics 
and  has  been  used  to  predict  the  combustion  behavior  for  axial  inlet, 
side-dump  inlet  and  bypass  configurations. 
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NOMENCLATURE 


A 


8 


C 


P 


E 


h 


constant  In  regression  rate  equation  (7) 
blowing  parameter 

specific  heat  at  constant  pressure 

activation  energy 
enthalpy 

heat  of  formation 

turbulence  kinetic  energy 
mass  flow  rate 
average  molecular  weight 
mass  fraction  of  species  i 

pressure 

heat  transfer  flux 
radial  distance 
fuel  regression  rate 
universal  gas  constant 
reaction  rate  of  species  i 

temperature 
axial  velocity 
radial  velocity 
tangential  velocity 

effective  transport  coefficient  of  species  i  =  u/o. 

incremental  distance  from  wall 
turbulence  dissipation  rate 
effective  viscosity 
density 

Prandtl  or  Schmidt  number  for  species  1 
mixture  fraction  (mass  fraction  of  total  fuel) 
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subscripts 
a  -  air 

c  -  chamber 

CO  -  carbon  monoxide 

CO^  -  carbon  dioxide 

CH  -  intermediate  hydrocarbon 
f  -  solid  fuel 

fu  -  unburned  fuel  vapor 

g  -  within  the  fuel  grain 

-  hydrogen 
H20  -  water 

'  oxygen 

p  -  near  wall  node  point,  port 

Tw  -  wall  temperature 

rad  -  radiation 

w  -  fuel  wall 
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I.  INTRODUCTION 


The  solid-fuel  ramjet  (SFRJ)  most  often  consists  of  a  solid-fuel  grain 
which  provides  the  walls  for  the  combustion  chamber.  A  sudden  expansion  at 
the  air  inlet  (either  axial  or  side-dump)  can  be  used  to  provide  flame 
stabilization.  Downstream  of  the  flow  reattachment  a  turbulent  boundary  layer 
develops  and  includes  a  diffusion-controlled  flame  between  the  fuel-rich  zone 
near  the  wall  and  the  oxygen-rich  central  core.  Due  to  that  diffusion  flame, 
heat  is  transferred  by  convection  and  radiation  to  the  solid  surface,  causing 
vaporization  of  the  fuel. 

At  the  Naval  Postgraduate  School  both  mathematical  modeling!-**  and 
experimental  efforts^- *0  have  been  conducted  to  determine  the  effects  of 
design  and  operational  variables  on  the  obtainable  performance. 

The  computer  model  simulation  of  the  SFRJ  combustion  process  has  evolved 
from  an  original  stream-function  vorticity  formulation*  to  a 
primiti ve-variable  (pressure,  velocity)  model  which  includes  an  aft  mixing 
chamber.3  These  axisymmetric  models  have  not  included  radiative  heat  transfer 
to  the  fuel  surface.  This  has  limited  the  utility  of  the  model,  since  many 
al 1 -hydrocarbon  fuels  produce  significant  amounts  of  radiative  transfer 
through  the  generation  of  carbon  particulates  in  the  flame  zone.  Recently4,  a 
simplified  radiative  model  was  included  in  the  primitive-variable  model  and 
resulted  in  improved  prediction  of  the  fuel  regressi on-rate  profile.  The 
primitive-variable  codes  were  based  on  the  Champion  2/E/FIX  computer  program 
developed  by  Pun  and  Spalding**.  The  Champion  code  is  a  two-dimensional  code 
and  therefore,  could  not  be  used  to  model  two  important  SFRJ  geometries; 
side-dump  and  bypass.  The  side-dump  configuration  has  air  injected  through 
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the  side  of  the  fuel  grain  and  the  dome  region  is  constructed  of  fuel.  The 
bypass  configuration  ducts  a  portion  (typically  501)  of  the  total  air  flow 
into  an  aft  mixing  chamber.  A  full  three-dimensional  (3-D)  model  is  needed  to 
predict  the  flow-field  in  these  configurations. 

The  side-dump  configuration  is  a  very  difficult  one  to  develop/optimize 
using  only  experimental  data.  There  is  a  complex  interaction  between  the 
location/size  of  the  inlet  dump,  the  length  of  the  dome  region  upstream  of  the 
dump,  and  the  resulting  fuel  regression  rate  pattern.  Regression  rate  data 
for  the  side  dump  geometry  have  not  appeared  in  the  open  literature.  However, 
several  general  observations  can  be  made.  Swirling  the  inlet  flow  can 
increase  and  smooth-out  the  downstream  fuel  regression  rate  pattern. 

Regression  rates  within  the  dome  region,  upstream  of  the  dump,  can  increase  or 
decrease,  depending  upon  the  dome  length,  number  of  inlets  and  amount  of 
swirl.  3-D  models  are,  in  this  case,  a  necessity  in  the  combustor  development 
process. 

Several  3-D,  elliptic,  primitive-variable  codes  have  been  developed.  One 
such  code,  developed  by  the  AiResearch  Manufacturing  Company, 12, 13  -jS  an 
extension  of  the  Champion  code.  The  former  computer  code  was  developed  to 
model  gas  turbine  combustors  and  formed  the  basis  of  the  3-D,  SFRJ  model 
discussed  herein. 
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II.  MODEL  OVERVIEW 


The  Garrett/AiResearch  model 
and  includes  finite-rate  chemical 
u»  v,  w,  h ,  k,  £,  <£ ,  m^ »  ^CH * 


is  3-D,  elliptic,  for  steady,  subsonic  flow 
kinetics.  The  dependent  variables  are 
and  m^g.  In  addition,  if  radiation  is 


included  (a  six-flux  model),  then  three  radiation  flux  vectors  are  calculated. 
Soot  emissions  can  be  calculated  if  desired.  Liquid  spray  calculations  can 
also  be  made  but  were  not  necessary  in  the  SFRJ  model.  The  turbulence  model 
is  the  same  as  used  in  the  Champion  program,  namely  the  modified  two-equation 
(k-e)  model  of  Jones,  Launder  and  Spalding^'^ 

Hydrocarbon  combustion  is  considered  to  be  a  four-step  process  as 
follows : 


CH  +  C  H  0  +  H„ 
x  y  x  y-2  2 


CxHy-2  +  \  °2 


xCO  +  ^2 


CO  +  1  o2  ♦  co2 

h2  +  j  o2  -  h2o 


(1) 

(2) 

(3) 

(4) 


Arrhenius  rate  expressions  are  used  for  each  of  the  above  reactions.  These 
chemical  kinetics  limited  rate  expressions  were  compared  to  rate  expressions 
based  on  an  eddy-breakup  model.  The  latter  expressions  attempt  to  account  for 
the  effects  of  turbulence  intensity  and  scale  and  species  concentration  on  the 
consumption  rate  of  a  particular  reactant.  The  two-part  boundary  layer 
(divided  by  y+  =  11.5)  used  in  the  Champion  code  is  also  employed.  Details  of 
the  computer  program,  the  differential  equations  and  the  line  by  line,  finite 
difference  solution  procedure  are  presented  in  references  12  and  13. 
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III.  MODEL  MODIFICATIONS  AND  SOLUTION  PROCEDURES 


The  major  modi f ications  required  of  the  Garrett/AiResearch  model  were 
geometric  and  the  inclusion  of  a  fuel  wall  on  which  finite-rate  chemical 
kinetics  could  occur.  The  geometric  changes  required  to  model  the  SFRJ 
configurations  were  readily  incorporated.  Several  additional  variables  and 
subroutines  were  necessary  to  include  the  regressing  fuel  surface  effects.  In 
addition,  fuel  property  subroutines  were  changed  to  reflect  the  properties  of 
Plexiglas  (PriM )  and  HTPB,  the  two  solid  fuels  used  in  the  present  study. 

The  governing  equations  on  the  boundary^  for  unity  Lewis  ,o er  are: 

Energy:  ^  =  (pv)w<W >  +  Fh  <!?>  w  (5) 

w 

(q  \  can  be  calculated  within  the  existing  model  m  was 
neglected  in  the  present  investigation.) 

3  m . 

Species:  R.  =  V.  (g-pl)  +  (pv)w(ni  -mi  )  (6) 

w  w  w  g 

(the  values  of  Rj  were  calculated  in  the  same  manner  as  for  the 
internal  nodes) 

It  was  assumed  that  the  fuel  regression  rate  could  be  represented  in  an 
Arrhenius  format:^-7 


pf?‘  =  (pv)w  =  A  pf  e'E/RTw  (7) 

Using  the  Couette  flow  approximation  for  the  boundary  layer  behavior  with  mass 
transfer!*^ 


*n(l+B) 
B=0  B 


(8) 


where  B  can  be  evaluated  from  the  enthalpy  as 


h  -h 

B  =  w 


h  -h-  -  (q  "  /fpv)  ) 
w  fg  'Vad  '  'w' 


(9) 


4 


(h)  calculate  (pv)  from  (5) 

(i )  calculate  Tw  from  (7) 

Iterate  (b)  to  (h)  until  (pv)w  converges 

( j )  calculate  pw  from 


R 

(k)  calculate  v,,  from  ^pv^w 

w  — 
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It  was  found,  as  previously  found  by  Netzer^  and  Stevenson  and  fletzer^, 
that  the  near-wall  turbulence  dissipation  had  to  be  increased  on  the  step  face 

(£p=kp^/^/0.46 )  and  that  the  grid  spacing  adjacent  to  the  fuel  surface  had  to 

be  maintained  slightly  less  than  that  required  for  y+  =11.5  in  order  to 

P 

obtain  temperature  distributions  in  qualitative  agrement  with  experimental 
data. 

In  order  to  obtain  convergence  it  was  necessary  to  input  initial 
conditions  which  were  in  qualitative  agreement  with  the  final  solution  (for 
example,  the  fuel  mass  fraction  had  to  be  specified  as  decreasing  from  the 
wall)  and  under  relaxation  (typically  0.1  to  0.5)  had  to  be  used. 

The  accuracy  of  the  solution  for  the  entire  combustor  of  the  solid  fuel 
ramjet  is  particularly  sensitive  to  the  boundary  conditions  on  the  fuel  wall. 
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IV.  RESULTS  AMO  DISCUSSION 


Four  geometries  were  investigated  vfig.  I);  axi symmetric  with  axial 
inlet,  bypass  with  axial  upstream  inlet  ,  single  side-dump  and  dual  side-dump. 
All  cases  included  an  aft  mixing  chamber.  For  the  axisymnetric,  axial -inlet 
geometry  the  predicted  reattachment  "points"  (where  u  near  the  wall  was 
interpolated  to  zero)  are  shown  in  Fig.  2.  The  reattachment  “points"  were  in 
good  agreement  with  the  maximum  heat  transfer  locations  measured  by  Krall  and 
Sparrow^®  but  shorter  than  those  measured  by  Phaneuf  and  Netzer®  in 
nonreacting  flow  with  wall  mass  addition. 

Normalized  regression  rate  profiles  are  shown  in  Fig.  3.  The 
experimental ly  measured  profiles  for  PMM  fuel  were  very  sensitive  to  inlet 
turbulence/distortion.  However,  the  3-D  code  (with  finite  rate  kinetics) 
predicted  results  very  nearly  identical  to  the  mixing  limited,  2-0  code 3, 
i.e.,  the  maximum  regression  rate  was  predicted  to  occur  upstream  of  the 
experimentally  measured  point.  The  average  regression  rate  was  in  good 
agreement  with  experiment. 

Examples  of  typical  radial  profiles  for  u,  nu  ,  mN  ,  and  T 

2  2 

near  the  grain  exit  are  shown  in  Fig.  4.  The  effects  of  including  finite-rate 
kinetics  and/or  the  injection  of  O2  that  occurs  into  the  boundary  layer  near 
the  reattachmnt  zone  are  readily  apparent;  a  broad  zone  of  high  temperature 
near  the  maximum  value  and  finite  concentrations  of  oxygen  below  the  flame. 

Figure  5  shows  regression  rate  profiles  at  three  axial  locations  for  the 
single  side-dump  geometry  without  inlet  air  swirl.  Mon-uniform,  high 
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regression  rates  occur  in  the  dump  region.  Lower  regression  rates  in  the  dome 
region  and  more  uniform  regression  rates  near  the  grain  exit  are  also 
observed. 

Figure  6  shows  the  predicted  regression  rate  results  with  two,  180° 
opposed  side-dumps,  with  and  without  inlet  air  swirl.  The  effects  of  inlet 
swirl  on  this  particular  geometry  are  readily  apparent.  Inlet  swirl  decreased 
the  regression  rates  in  the  dome,  but  significantly  increased  the  rates 
downstream.  Regression  rate  uniformity  was  also  significantly  improved  except 
adjacent  to  the  inlet  dump. 

Figure  7  shows  the  results  of  increasing  the  dome  length  and/or  the  inlet 
dump  area.  Increasing  the  inlet  area  resulted  in  increased  fuel  consumption 
in  the  dome  region  and  increased  nonuniformity  downstream.  Lengthening  the 
dome  above  the  reference  value  signi f icantly  reduced  head-end  fuel  regression 
rates  without  major  effects  in  the  downstream  regions. 


V.  CONCLUSIONS 


The  model  predictions  show  that  significant  variatiions  in  fuel 
regression  rate  (both  circumferentially  and  axially)  are  caused  by  the  manner 
in  which  the  air  is  injected  into  the  fuel  grain.  In  general,  inlet  swirl 
increases  and  smooths-out  the  fuel  regression  rate  downstream  of  a  side-dump 
inlet.  The  effects  of  inlet  location  and  size  on  the  regression  rate  behavior 
are  difficult  to  estimate  a-priori  for  a  specific  set  of  test  conditions.  The 
model  must  be  used  to  predict  the  expected  results.  Much  additional  work  is 
required  in  order  to  validate  the  predictions  of  the  3-0  SFRJ  model.  However, 
the  ability  to  examine  the  effects  of  inlet  air  location  and  flow 
characteristics  on  the  fuel  regression  rate  pattern  (and  the  combustion 
efficiency,  etc.)  has  provided  a  needed  improvement  to  the  earlier  2-0  model. 

Current  work  is  being  directed  at  incorporation  of  the  soot 
production/consumption  and  radiation  subroutines  into  the  SFRJ  configuration 
and  combustion  environment.  Regression  rate  profiles  for  the  more  recent 
configurations  are  also  being  compared  to  experimental  data  as  they  become 
available. 
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Figure  1.  Solid  Fuel  Ramjet  Geometries  in  Model. 


Figure  2.  Reattachment  Locations  for  Axial-Inlet  Flow. 
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Figure  3.  Fuel  Regression  Rate  Profiles  for  Axial-Inlet  Flow 
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Figure  4.  Predicted  Radial  Profiles  for  Gas  Properties  Near  the  Grain 
Exit,  Axial -Inlet  Flow. 
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Figure  5.  Fuel  Regression-Rate  Profiles,  Single-Side  Dump,  No  Inlet  Swirl. 


Figure  6.  Effects  of  Inlet  Swirl  on  Fuel  Regression  Rate  Profiles,  180°  Opposed 
Side- Dumps. 
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Opposed  Side-Dumps  with  Inlet  Swirl. 
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